Constrained dynamics of a polymer ring enclosing a constant area 
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The dynamics of a polymer ring enclosing a constant algebraic area is studied. The constraint of 
a constant area is found to couple the dynamics of the two Cartesian components of the position 
vector of the polymer ring through the Lagrange multiplier function which is time dependent. The 
time dependence of the Lagrange multiplier is evaluated in a closed form both at short and long 
times. At long times, the time dependence is weak, and is mainly governed by the inverse of the 
first mode of the area. The presence of the constraint changes the nature of the relaxation of the 
internal modes. The time correlation of the position vectors of the ring is found to be dominated by 
the first Rouse mode which does not relax even at very long times. The mean square displacement 
of the radius vector is found to be diffusive, which is associated with the rotational diffusion of the 
ring. 

PACS numbers: 36.20.-r Macromolecules and polymer molecules - 05. 40. -a Fluctuation phenomena, random 
processes, noise and Brownian motion - 02.50.Ey Stochastic processes - 05.40.Jc Brownian motion 



I. INTRODUCTION 

The equilibrium statistics of a planar Brownian mo- 
tion enclosing a constant area was first studied by Levy 
in 1940 with a view to obtain theprobability distribution 
function for the algebraic area 0, 0| ■ A simple analogy 
to topologically constrained (entangled) polymers, which 
can approximately be represented by a planar random 
walk constrained to enclose a constant area, was later 
pointed out by Brereton and Butler §j. To describe 
the configurational and mechanical properties of these 
polymers, they re-addressed this problem using discrete 
Gaussian chains. An exact expression for the probabil- 
ity distribution function was subsequently obtained using 
continuous models for Gaussian chains jj, Q . 

Most of the studies on constrained polymer rings have 
been restricted to static equilibrium behavior. However, 
the dynamics of a polymer ring constrained to enclose 
a constant area can serve as a useful model in several 
different areas of statistical physics. The Rouse dynam- 
ics of a highly dense ensemble of long polymer chains, in 
general, is a simple and well known example of the con- 
strained dynamics due to the presence of entanglements. 
The highly non-local nature of these constraints renders 
the analytical solution of the problem difficult. There- 
fore, the dynamics of polymer melts is mainly studied by 
phenomenological tube models |(| and their extensions 

7]. For a polymer ring, on the other hand, analogies 
with topologically constrained loops can be realized at a 
more formal level 0. The mathematical formulation of 
the constant area constraint is strongly related, for in- 
stance, to the Gaussian invariance of the entangled loops 

6] , and topological problems concerning the writhe of a 
polymer ring |8(. Some connections with rings enclos- 
ing topological obstacles can also be discussed within a 
similar framework 0]. 

As another instance, the polymer ring can be regarded 
as a very simple model to describe the dynamics of a two- 
dimensional vesicle bounded by a one-dimensional lipid 



membrane, the perimeter of which undergoes randomly 
kicked motion due to the stochastic thermal forces from 
the solvent; the osmotic pressure difference between the 
two sides of the membrane is such that it keeps the area 
constant at all times. The constraint of constant area 
can be enforced by introducing a time dependent La- 
grange multiplier in the stochastic equations for a poly- 
mer ring. The Lagrange multiplier, therefore, plays the 
role of the time dependent differential pressure, the dy- 
namics of which ensures that the area is constant at all 
times. 

The technical difficulties involved in treating the con- 
strained dynamics of a polymer ring arise from the global 
nature of the area constraint. Additionally, since the con- 
straint of a constant area has to be satisfied at each time 
t, the dynamical behavior of the two Cartesian compo- 
nents of a spatial vector r(s, t) = (r x (s,t),r y (s,t)) at any 
point s along the curve are no longer independent — the 
rigid constraint of a constant area strongly couples the 
two coordinates through the Lagrange multiplier. The 
time dependence of the latter makes the analytical solu- 
tion of the stochastic equations difficult. 

In a more general context, the present problem belongs 
to the comparatively less understood field of stochastic 
dynamics with rigid constraints. Even in a more specific 
case of the dynamics of polymers with rigid constraints, 
the literature is mainly restricted to the Kirkwood theory 
for hydrodynamics and the dynamics of semiflcxible poly- 
mers One reason for the limited literature on such 
systems is the technical difficulty involved in dealing with 
dynamical problems with rigid constraints. The nature of 
the difficulty become apparent in the field-theoretic for- 
mulation of the dynamical problem, where the presence 
of any constraint changes the number of dynamical vari- 
ables and adds the Lagrange multiplier to the set of vari- 
ables. Any functional formulation of the dynamics of the 
constrained system is, therefore, altered by an additional 
Jacobian, which is difficult to determine. In some cases 
the introduction of non-commuting Grassmann variables 
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can facilitate a solution ^(j- F° r the present problem, 
however, we show that the nature of the area constraint 
is such that a simple Langevin description suffices and a 
field-theoretic formulation is not neccesary. 

Within the Langevin description, the area constraint is 
accounted for by the Lagrange multiplier function which 
is time dependent. The time dependence of the Lagrange 
multiplier can be calculated in a closed form (which is 
quite unusual for the problems of this nature) for an ini- 
tial perturbation of the area modes and allows us to solve 
the coupled stochastic equations for the polymer ring. It 
is shown that the dynamics of a polymer ring enclosing 
constant area is mainly governed by the lowest Rouse 
mode. The constraint of a constant area is fixed by the 
first Rouse mode (associated with rotation 6]), which is 
found to be diffusive at long times. The present problem 
is a simple and well chosen example to describe several 
interesting aspects of stochastic dynamics with rigid con- 
straints. 



II. MODEL 

Let us first introduce the model we are going to study 
and review some of the static results. The algebraic area 
enclosed by a planar random walk in a continuous repre- 
sentation is defined by 



A{ r (s,t)} = ^J da Lt)x^| k. 



ds 



(1) 



where k is the unit vector perpendicular to the (x, y) 
plane; r(s, t) is the position vector starting from the cen- 
tre of mass of the chain (placed at the origin) to the 
segment s along a chain of contour length N , which for 
a ring satisfies r(N, t) = r(0, t). In statics, the constraint 
of a constant area is imposed by including a delta func- 
tion <5(A{r} — A) in the partition function such that the 
probability distribution function for A is given by 
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ds 



dr(s) 
~ds~ 



, (2) 



where TV is the normalisation constant and b is the seg- 
ment length; the term in the exponential is the Wiener 
measure which accounts for the entropic elasticity of the 
polymer. The area constraint restricts the conformations 
of a polymer ring in such a way that closed random walks 
of only a given constant area A contribute to the parti- 
tion function. Since the sign of the area (as defined by 
Eq. (JIJ) depends on the orientation of the contour, the 
constraint restricts the magnitude of the area but has no 
control over the shape of a ring. The latter implies that 
a random walk with a large number of wiggles along its 
contour length will undergo many cancellations to keep 
the average area constant 3]. 



The delta function in Eq. can be Fourier trans- 
formed, i.e. S(A{r} - A) = J dg exp(ig(A{r} - A)), 
such that g is a variable conjugate to A. The prob- 
ability distribution P(A, N) can then be determined 
for a discrete closed random walk using normal mode 
analysis [3j or its continuous representation using func- 
tional integration 0, HJ; it is given by P(A, N) — 
(2Nb 2 cosh 2 (2nA/ Nb 2 ))^ 1 , which has the asymptotic be- 
havior of exp(-47rA/M> 2 ) for A/Nb 2 > 1. 

The constrained dynamics of a constant area can be 
treated by the method of Lagrange multipliers @; the 
idea is to add the constraining force explicitly in the 
Langevin equation. The constraint of a constant area 
when introduced using a Lagrange multiplier function 
X(t) produces a coupling between x and y coordinates. 
The coupled Langevin equations are given by 
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dr x (s,t) 
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> dr v (s,t) 

' at 



2k B T d 2 r x (s,t) ^ M dr y (a,t) 
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where the first term in the above equations represents 
dissipation due to the viscous forces of the solvent with 
friction coefficient (; f is the fluctuating force whose 
time average is zero. The dissipation and the fluctu- 
ations are related to each other by the simplest form 
of the fluctuation dissipation theorem which dictates 
(/ l (s,t)/ J (s',i / )) = 2Ck B T8 ij §{s-s')8{t-t l ). Thesecond 
term represents the elastic force due to the chain connec- 
tivity. The coupling shows the most trivial effect of the 
constraint: the dynamics of the two components r x {s,t) 
and r y (s, t) is no longer independent, as is the case in the 
absence of constraint. 

In terms of the normal modes, i.e. r p(t) = 

i f*ds r(s, i)e 2ipTrs/JV , Eqs.Q and gj can be rewrit- 
ten as 



dr px (i) 


-8p 2 ir 2 k B T 


dt 


Nb 2 




+ fpx(t) 


dr py (t) 


-8 P 2 TT 2 k B T 


dt 


Nb 2 




+ /py(0i 



r px (t) - 2ipnk B T\(t)r py (t) 

(5) 

r py {t) + 2ipTrk B TX(t)r px (t) 

(6) 



where Co = N( and C P #o = 2AT£6]. A(t) = 
E P l^oo(-^ 7r ) yl p(*)' where 

A p (t) = rpx(t)r-py(t) - r py (t)r- px (t). (7) 

A more compact vectorial representation of the set of 
coupled Langevin equations can be presented by using a 
matrix notation, 



<p|r p (t) 



-k B TM p (t)-r p (t)+f p (t) 



(8) 
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where the matrix M p (i) is defined 



M„(t) = 



2 n 2 /(Nb 2 ) 2ipirX(t) 



-2ipnX(t) 



2 7T 2 /(Nb 2 



whose eigenvalues £ e are given by 

4 = 8p 2 ir 2 /(Nb 2 ) ± 2p7rA(t), 



(9) 



(10) 



which correspond to the relaxation of the individual 
Rouse modes of the chain the details of which will be 
discussed later. 

Since the area constraint is local in time t but global in 
the contour length s, the Lagrange multiplier A depends 
only on time. The time dependence of Lagrange multi- 
plier makes the exact solution of the coupled Langevin 
equations difficult. Any exact or approximate solution 
requires a complete knowledge of the time dependence of 
A(t). Since the dynamics of A(i) has to be such that the 
area A(t) has to be constant at all times, we begin by 
writing down the equation of motion for the area modes 
itself, the details of which are presented below. 



III. AREA CONSERVATION 

A. Dynamics of mode-dependent area 

To get a dynamical equation for the area of a poly- 
mer ring, we add and subtract the imaginary part of Eq. 
© from Eq. JSJ in such a way that we get dynamical 



equations for (r p 



t ) and (r_ 



A series of simple steps will provide a stochastic 
(Langevin type) equation for area: multiply the dynami- 



cal equation for (r 



px 



' I'll 



) by (r_ 



px 



-I'll 



) (that is its 



complex conjugate) and (r- px - ir_ py ) by (r px + ir py ); 
after having added the two equations, separate the real 
and the imaginary parts of the resulting equation, which 
for p =/= modes are given by the following expressions 
respectively: 



Cp 



dr p 2 {t) 

dt 
dA p (t) 

dt 



= - {2kp 2 + Apnk B TX{t))r p 2 {t) (11) 
= -(2kp 2 + 4 P TTk B T\(t)) A p (t), (12) 



where k = 8ir 2 k B T/Nb 2 , A p (t) is given by Eq. Q and 
r p 2 (t) = r px (t)r- px (t) + r py (t)r- py (t). It is not surpris- 
ing that the above two equations have a similar structure 
since the first one represents the radius and the second 
one the area of the pth mode of the ring. To get the 
time dependence of A(<), we are interested in the time 
development of the average of A p (t), which in thermal 
equilibrium is time independent. Therefore, in writing 
the dynamical equation for A p (t), we have implicitly as- 
sumed an average over a non-equilibrium ensemble in 
which (A p (t)} is explicitly time dependent. This gives 
us the dynamical response of the Lagrange multiplier to a 



perturbation of the area modes away from their thermal 
equilibrium values |ll|. 

The area conservation demands dA(t)/dt = 
^2 p (—ipir)dA p (t)/dt = 0. This constraint when 
used in Eq. (|12fl produces the following expression for 
the Lagrange multiplier: 



A(t) 



4tt E P p 3 Mt) 

Nb 2 Y. P P 2 MtY 



(13) 



where Nb 2 /Att is the area of the first mode; it corresponds 
to a random walk which on an average forms a circle. 
For a random walk to form a circle in an average sense, 
the circumference 2nr should be equal to N 1 / 2 b, where 
r is the radius of the circle; this condition determines 
r = N 1 / 2 b/2ir. Thus, the area of the circle on an average 
turns out to be Nb 2 /An. The second mode is a symmetric 
figure of eight and so on. 

The time dependence of the Lagrange multiplier can 
be determined in a self-consistent way. To begin with 
we consider X(t) to be independent of time and given by 
— Air/Nb 2 . This approximation amounts to solving Eq. 
(|12fl . which is a simple first order differential equation 
in time, and can be easily solved to produce A p (t) — 
A p (0) exp(— t/r p ), where t p is the relaxation time of the 
pth mode given by t p = N 2 b 2 (/[8TT 2 k B Tp(p - 1)]. The 
latter implies that all modes except the first one, which 
has a constant area Ai(0), decay with the relaxation time 
of t p . The same is true for the radius of the pth mode. 

We can now substitute the expression for A p {t) into Eq. 
(|13fl to see how X(t) evolves with time. In the limiting 
case of X(t — > oo ) = —Air/Mb 2 , the dominance of the first 
mode is evident. The latter is true for all times greater 
than the relaxation time r = N 2 b 2 ^/[8ii 2 k B T]. It is clear 
from Eq. (|13fl that for all times t <C r, the dynamics 
is determined by the summation over large number of 
modes; for times t 3> r, on the other hand, the dynamics 
is dominated by the first mode. Thus, the dynamics of 
X(t) is reflected in its dependence on p modes, which 
decays in time in such a way that for t 3> r it is only 
determined by p = 1. 

To get an idea of the time dependence for times much 
less than the relaxation time r, one can replace the sum- 
mation by an integration over p modes in Eq. l|12fl . 
In the limit of ( « t, the time dependence of the 
Lagrange multiplier is approximately given by X(t) w 
—AiT/Nb 2 (T/t) l l 2 . It is to be noted that in carrying out 
the integration over the modes, we have neglected the p 
dependence of A p (0). The mode dependence of A p (Q), 
however, does not alter the way X(t) scales with t. At 
short times, the nature of the decay of the pth mode 
is reflected in the time dependence of X(t). It is evi- 
dent from the expression for A p (t) that for t r/p 2 , 
the pth mode has already relaxed and does not con- 
tribute to the summation in Eq. (|13fl . This implies 
that the number of modes that contribute to the sum- 
mation in Eq. i|13|) at any time t are approximately given 
by P*(t) ~ (r/i) 1 / 2 . The algebraic time dependence of 
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X(t) can easily be understood from Eq. (|12fl . If it is as- 
sumed that in this time regime some modes have already 
relaxed, then dA p * jdl = for these modes. Then the so- 
lution for the Langrange multiplier X(t) is approximately 
given by 



X(t) 



47T 
iVfc 2 



V 



(14) 



Since the modes follow the Rouse relaxation, the inverse 
square root time dependence of X(t) can be understood. 

The time dependent expression for X(t) when substi- 
tuted into Eq. (O yields A p {t) « A p (0)e-P 2 */ T eP(*/ T ) 1/2 , 
which at short times is a stretched exponential. The lat- 
ter implies that for time t C t, the constraint of con- 
stant area affects the dynamics of A p {t) in such a way 
that the lower modes grow at the expense of the higher 
ones leaving the total area constant; for time t r only 
the first mode, which attains a constant value, survives. 
In contrast, the dynamics of X(t) is such that it decays 
in time as (r/t) 1 / 2 and at t ^> r attains a constant value 
of -iir/Nb 2 . 

This gives us an idea about the dynamics of A(i), which 
for all times t 3> r, is weakly dependent on time and is 
mainly dominated by the first mode. In what follows, we 
will replace X(t) by a time independent quantity Ao to 
discuss the statics of a closed random walk. 



B. Statics 

To recover some of the equilibrium results of the earlier 
studies from dynamics, we begin by introducing temporal 
Fourier transform into Eqs. (JSJ and ©: 

- iuiC, p r px (uj) = -kp 2 r px (uj) - 2ipnX k B Tr py (uj) 

+ JWM (is) 

-iw( p r py (w) = -kp 2 r py {ui) + 2ipnX k B Tr px (uj) 

+ /«/("). (16) 

The resulting coupled equation can be easily be solved 
to calculate the following correlation: 

(r p (w)T_ p (-w)) = 2k B TQ p 



(C 2 Lo 2 + (kp 2 -2pirk B TX a ) 2 ) 
1 

{Quj 2 + (kp 2 + 2pwk B TX ) 2 )\ ' 

(17) 



The mean square distance of a chord starting from 
r(0) and ending at any point s along the chain con- 
tour is defined by R 2 (s) = (r(s) — r(0)) 2 . In terms of 
the normal modes it is given by R 2 (s) = 4^^ 1 (r p • 

r_ p ) sin 2 (p7rs/iV). In the limit of N — > oo andpn/N — » q, 
the summation can be converted to an integral to give 



R 2 ( s ) 



r° 

2b 2 /tt / 
Jo 



dq S m 2 (qs)/(q 2 + g 2 b A /16), (19) 



which is the same as Eq. (5.3) in the paper by Brereton 
and Butler 3J. In writing the above expression we have 
taken Ao = ig such that the integral can be done by 
contour integration; it has double poles in the complex 
g plane given by g — ±4iq/b 2 |3J. The latter implies 
that when q = pir/N, the Lagrange multiplier is given by 
A = ±4: P n/Nb 2 . 

The integral in Eq. I|19|) can easily be evaluated and 
Fourier transformed with respect to <?; the result can 
be averaged over the equilibrium probability distribution 
given by P(A) = exp(— inA/Nb 2 ) to give the following 
expression: 



<R 2 W>, 



^ {1 + ^J^±fl) ,_„■, 



where a — 2A/Nb 2 . It is interesting to note that 
in equilibrium the probability distribution, which yields 
(A) = Nb 2 /4n, imposes the dominance of the first mode 
over all other modes. In essence it means that the pres- 
ence of the area constraint introduces significant devia- 
tions from the unperturbed random walk for A/Nb 2 3> 1; 
in the latter limit the mean square distance is given by 
(R 2 ( S )) A ^sb 2 (l + S /N) 2 /2. 



C. Dynamics 

As discussed in Section III, the dynamics of the La- 
grange multiplier function X(t), when calculated in a self- 
consistent way, suggests its inverse dependence on the 
first mode of the area, i.e., -4n/Nb 2 . The latter expres- 
sion when substituted into Eqs. I|15f) and jTSjl produces 
coupled equations in r p (x) and r p {y), which can easily be 
solved to give 



(r p (w)T_ p (-a;)> = 2fc B TC P 



1 



(k 2 p 2 (p + l) 2 ) 
1 



{Quo 2 + k 2 p 2 (p- l) 2 ) 



(21) 



The time-correlation function can be calculated using 



(r p (t) ■ r_ p (0)> = / % (r» • r_ p (- W )> e***. For p = 1 



where the angular brackets represent an average 
with respect to the random noise. To calcu- 
late the equal time correlation we use {r pl r- po ) = mode the correlation function is given by 
f (du>/2ir) (r p i(u>)r^ p j(—uj)); the resulting expression is 
given by 



Nb 2 



(r p (t) ■ r_ p (t)> = 



iV6 2 /2 



( P 2 tt 2 - AgiV 2 6 4 /16)' 



(18) 



16tt 2 
k R T\t 



exp 



knT\t - t'\ 



n Kb 



N 2 b 2 C 



(22) 
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For p ^ 1 modes, on the other hand, the time-correlation 
function is given by 



(r p (i)T_ p (i')) 



Nb 2 ( An 2 p 2 k B T\t-t'\ 

exp ' 



4tt 2 (p 2 - 1) * V A^C 
, ( Air 2 pk B T\t - t' 

C0Sh { N 2 b 2 C 

1 . , ( A-K 2 pk B T\t-t'\ 
H — sinn 1 
P 



V A 2 6 2 C 



(23) 



The radius of the ring at any point s along the curve is 
given by R(s, t) — r(s,i) — i"o(i), where ro(t) is the po- 
sition of the centre of mass at time t defined as ro(t) = 
f Q dsr(s,t). The time correlation of the radii at two 
different times in terms of the normal modes is given by 

(R(M)-R(M')) = 4££Li(r P (*)T-p(f)>- Since the 
relaxation time associated with the first mode, that is 
t = N 2 b 2 (/8Tr 2 k B T, is the longest relaxation time of 
the correlation function, it corresponds to the rotational 
relaxation time @- For all times greater than the re- 
laxation time, that is t 3> t, the dynamics is dominated 
by the first Rouse mode, which is diffusive. The mean 
square displacement of the radius R( A, t) is given by 



((R(A,t)-R(A,0)) 2 ) = 



8k B Tt 



(24) 



where at long times it corresponds to the rotational diffu- 
sion of the ring; the rotational diffusion constant is given 
by D rot — 2k B T/N(. The mean square displacement 
of the centre of mass, on the other hand, is estimated 
by the p = mode, i.e. ((r (i) - r (0)) 2 ) = ^P; 
at long times the self diffusion constant of the centre 
of mass of the ring is given by D Q = k B T/NQ 6]. The 
rotational diffusion constant is, therefore, twice of the 
diffusion constant of the centre of mass of the ring. It is 
not surprising, however, that both the diffusion constants 
have Rouse type scaling since the only way the diffusion 
constant can depend on the system properties is through 
D oc k B T/((b 2 /A)®. The area is given by A cx (l/b 2 N); 
the diffusion constants for the p = diffusion and the 
area preserving p = 1 rotation must, therefore, have the 
same Rouse type scaling. 

Since the area constraint couples the x and y modes, 
the cross correlations, i.e. (r px (t)r-py(t')), turn out to be 
non-zero. The equal time correlation amounts to estimat- 
ing the average of A p (t), which can easily be evaluated 
using Eq. (J7J . The summation over p modes defines the 
average area, which is given by (A(t)) = Nb 2 /An. 

The dominance of the diffusive first mode has been seen 
in a very different context of grafted polymer brushes 
[l2| . In the latter study the distribution of the chain ends 
is accounted for by self-consistent-field theory which de- 
termines the mean-field potential seen by each monomer 



self-consistently |l3§ . In dynamics it amounts to the dom- 
inance of the first mode, which does not relax, and is as- 
sociated with the diffusion of the chain end in brushes. 
In the present study although the imposition of area con- 
straint through the time dependent Lagrange multiplier 
and its determination from the dynamical equations is of 
very different nature, the dominance of the first (Rouse) 
mode, as also seen in the dynamics of polymer brushes, 
is interesting. 



IV. CONCLUSIONS 

In summary, the present theory illustrates the effects 
of a relatively simple global constraint on the stochastic 
dynamics of non interacting polymer chains. We show 
that the presence of the area constraint couples the dy- 
namics of the Cartesian components of the spatial vector 
r(s,i) = (r x (s,t),r y (s,t)) through the time dependent 
Lagrange multiplier. The time dependence of Lagrange 
multiplier can be determined as a dynamical response 
to the initial perturbation of area modes such that the 
condition dA(t)/dt = is satisfied at all times. At very 
short times, that is t <C t, the Lagrange multiplier de- 
cays slowly with an inverse dependence on time, where 
r = N 2 b 2 (/8TT 2 k B T. For long times, that is t > t, the 
Lagrange multiplier is mainly determined by the inverse 
of the first mode of the area, which on an average is a 
random walk about a circle. As opposed to the relax- 
ation of the internal modes of a ring without constraint, 
the area constraint changes the nature of relaxation of 
these modes dramatically. At long times the dynamics is 
completely determined by the first (Rouse) mode which 
does not relax. The first mode is associated with the ro- 
tation of the ring; the mean-square displacement of the 
radius at long times is found to be diffusive. 

Although the idea was to present a general theory for 
the dynamics of area preserving systems, the present 
approach can serve as a simple description to under- 
stand some of the dynamical aspects of polymer sys- 
tems under topological constraints, entanglements and 
two-dimensional vesicles — systems where the area con- 
servation has to be explicitly accounted for. 
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